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Abstract: If a one-dimensional quantum lattice system is subject to one step 
of a reversible discrete-time dynamics, it is intuitive that as much "quantum 
information" as moves into any given block of cells from the left, has to exit 
PL|' that block to the right. For two types of such systems — namely quantum walks 

and cellular automata — we make this intuition precise by defining an index, 
a quantity that measures the "net flow of quantum information" through the 
system. The index supplies a complete characterization of two properties of the 
discrete dynamics. First, two systems Si^S? can be "pieced together", in the 
sense that there is a system S which acts like S\ in one region and like S2 
in some other region, if and only if S\ and S2 have the same index. Second, 
' the index labels connected components of such systems: equality of the index is 

necessary and sufficient for the existence of a continuous deformation of Si into 
S2. In the case of quantum walks, the index is integer- valued, whereas for cellular 
automata, it takes values in the group of positive rationals. In both cases, the 
map S i-> ind S is a group homomorphism if composition of the discrete dynamics 
is taken as the group law of the quantum systems. Systems with trivial index are 
precisely those which can be realized by partitioned unitaries, and the prototypes 
of systems with non-trivial index are shifts. 



1. Introduction 



Quantum walks and quantum cellular automata are quantum lattice systems 
with a discrete step dynamics, which is reversible, and satisfies a causality con- 
straint: In each step only finitely many neighboring cells contribute to the state 
change of a given cell. This leads to an interesting interplay between the condi- 
tions of reversibility (unitarity) and causality, which is the subject of this article. 

Starting point of the analysis is a simple intuition: for any connected group 
of cells in a one dimensional system as much "quantum information" as moves 
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into the subsystem from the left has to move out at the other end. Moreover, 
this "flow" is a conserved quantity, in the sense that it remains constant over 
the spatial extent of the system. It can thus be determined locally at any point. 

Making this intuition precise, we associate with every such lattice system an 
index, a quantity measuring the net flow of information. The index theory devel- 
oped in this work completely resolves three, a priori very different, classification 
problems: 

(1) Find all locally computable invariants. It is shown that there exists a 
"crossover" between two systems Si , S2 if and only if their indices coincide. More 
precisely, a crossover S between S\ and S2 is a system which acts like S± on a 
negative half line {x\x < a} and like S2 on a positive half line {x\x > b}. Clearly, 
a locally computable invariant must assign the same value to two systems if there 
exists a crossover between them. It follows that any invariant is a function of 
the index. 

(2) Classify dynamics up to composition with local unitaries. A natural way 
of constructing dynamics which respect both reversibility and causality is by 
concatenating layers of block unitaries. In every step, one would decompose the 
lattice into non-overlapping finite blocks and implement a unitary operation 
within every block. Such local unitary implementations are conceptually related 
to the gate model of quantum information. Not every time evolution may be 
realized this way: a uniform right-shift of cells serves as the paradigmatic counter- 
example. We show that the systems with local implementations are precisely 
those with trivial index. Consequently, equivalence classes of dynamics modulo 
composition with block unitaries are labeled by their indices. 

(3) Determine the homotopy classes. It is proved that two systems may be 
continuously deformed into each other (with a uniform bound on the causality 
properties along the connecting path) if and only if they have the same index. 

We will consider the above questions, and define indices, for two kinds of 
systems. Quantum walks are, on the one hand, the quantum analogs of classical 
random walks. On the other hand, they are discrete time analogs of a standard 
quantum particle "hopping" on a lattice according to a Hamiltonian which is 
a lattice version of the momentum operator \d x . The index defined for these 
systems is the same as a quantity called "flow" by Kitaev [T^]. Intuitively, this 
measures the mean speed of a quantum walk, expressed in units of "state space 
dimensions shifted to the right per time step" . The mathematical background 
has been explored, in a more abstract setting, by Avron, Seiler and Simon [3J. 
Kitaev's work treated the first classification problem above. We will re-prove 
his results with an eye on generalizations to cellular automata, and will supply 
solutions to questions (2) and (3). Although the quantum walks seemed to be 
comparatively straightforward initially, the intuition gained from this case served 
us well in setting up the theory for the much more involved case of cellular 
automata. This allowed us to build an abstract index theory covering both cases 
with almost identical arguments (Sects. [3] and [5]) . 

Cellular automata are characterized by the property that whatever state is 
possible in one cell (e.g., a superposition of empty/occupied) can be chosen in- 
dependently for each cell. Expressed in terms of particles this means that we 
are necessarily looking at a "gas" system of possibly infinitely many particles. 
The basic definition of quantum cellular automata was given in |I6j . On the one 
hand, the setting considered here is more restrictive than |16j , covering only one- 



Index theory of one dimensional quantum walks and cellular automata 



3 



dimensional systems. On the other hand, we are allowing for non-translationally 
invariant dynamics — a strong generalization over the earlier paper. In fact, 
having completed the present work we feel that the translation invariance as- 
sumed in |16j was obscuring the fundamental interplay between reversibility and 
causality. Accordingly, we obtain here a stronger structure result, even though 
it is built on the same key ideas. Throughout, there is a strong interplay be- 
tween local and global properties. For example, the following statement is an 
immediate consequence of our main Theorem [5] If a nearest neighbor cellular 
automaton has somewhere a cell of dimension n, and somewhere else a cell of 
dimension to, coprime to n, then it can be globally implemented as a product of 
two partitioned unitary operations. 

Our paper is organized as follows: After giving two examples in Sect. [31 We be- 
gin by a mathematical description of what we mean by quantum walks (Sect. [33} 
and quantum cellular automata (Sect. [32]). We then describe the notion of locally 
computable invariants, and why they should form an abelian group (Sect. [5]). 
A similar general explanation of the notion of local implementation is given in 
ScctU] The detailed theory for quantum walks is in Sect. [HI and in Sect. [7] for 
cellular automata. This includes the proof that an index previously defined in 
the classical translationally invariant case [10] . coincides with our index for this 
special case (Sect. [73)l . We close with an outlook on variants of index theory for 
either higher dimensional systems or automata with only approximate causality 
properties (Sect. [5}. 

2. Examples 

Before introducing the mathematical setting, we will illustrate the problems 
treated in this paper by giving two concrete examples. 

2.1. Particle hopping on a ring. The simplest example is given by a single par- 
ticle on a ring of N sites arranged in a circle. More precisely, the Hilbert space 
we are considering is C N with basis vectors {|eo), ■ • ■ , |ejv-i)}. The vector \a) is 
taken to represent a "particle localized at position i" . One step of a reversible 
discrete-time dynamics is simply given by an arbitrary unitary U € U(C N ). 
We will consider two such time evolutions Uo,Ui defined with respect to the 
standard basis by 

) (1) 

respectively. The first unitary is the trivial evolution and the second one models 
a uniform movement of the particle with velocity one site per time step to the 
right. 

The physical interpretation of this simple model seems clear: we can think 
of a lattice version of a particle with one spatial degree of freedom, where we 
have introduced cyclic boundary conditions to get a simple, finite description. 
The causality property defining a quantum walk then expresses the physically 
reasonable assumptions that couplings are local and dynamics preserve locality. 

In this setting, it is natural to think of the time evolution as being generated 
by a Hamiltonian: Ut = e ltH . Such a Hamiltonian formula would allow us to 
extend the dynamics to arbitrary real times t € 1Z. To recover H, we need to 
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Fig. 1. Absolute value squared of a vector initially localized at site 6 (of 10) under the action 
of the time evolution e ltH for t g {0, .25, .5., .75, 1}. The Hamiltonian is given by Eq. J2J. The 
dynamics is not causal for non-integer values of t. 

take a logarithm of U\. This operation is of course not uniquely defined, but the 
ansatz 

ff=£^i/*X/*l ( 2 ) 

in terms of the Fourier basis 

AT-l 
3=0 

seems particularly appealing. It can easily be checked to be compatible with our 
previous definitions of Uq and U\. 

Docs this definition yield physically satisfactory dynamics for all times t £ 
[0, 1]? Hardly. As demonstrated in Fig. (TJ a particle initially localized at site 6 
(N = 10) will spread out over the entire ring during the interval t £ [0, and 
rcfocus to site 7 during t £ [i, 1]- Any sensible notion of "causality" is violated for 
non-integer values of t. Note that this contrasts with the time evolution generated 
by the momentum operator H = id x of a continuous- variable system. The latter 
fulfills (e ltH ip)(x) = ip(x — t), thus preserving the localization properties of 
vectors tp £ L 2 (R) for any t £ TZ. 

The discussion immediately raises several questions: 

Is there a way to continuously interpolate between Uq and U\ while preserving 
causality? Our particular choice © for H could have been unfortunate. Con- 
ceivably, there is a better-suited, potentially time-dependent Hamiltonian which 
docs the job. More physically phrased: is it possible to discretize the spacial 
degrees of freedom, but not the temporal ones, of a free single particle, while 
respecting causality? (The answer is: no, this is not possible). 

Is there a simple way of deciding that the answer to the first question is 
negative? (Yes: the index of Uq is 0, the index of ~U\ is 1. The index is constant 
on connected components with uniformly bounded interaction length). 

Further questions we will answer include: Is there a "gate-model" implemen- 
tation of Z7i? More precisely, can I write U\ in terms of a constant-depth sequence 
of unitary operations, each of which acts non-trivially only on a constant num- 
ber of sites? (No). Can I engineer a large system on JV' » f sites, endowed 
with a global time evolution U, such that the restriction of U to some region of 
contiguous sites looks like Uq and the restriction of U to some other region looks 
like U x l (No). 

All these statements are made precise in Theorem [3] 
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2.2. Cluster state preparation. In this section, wc consider ./V spin-1/2 particles 
arranged on a circle. The Hilbert space of the system is H = ® !=1 C 2 . In the 
previous example, it was clear what it meant for a state vector ip to be "localized" 
in some region A C [1, N]: namely this was the case if (t/j\si) = for all i £ A. 
For state vectors on tensor product spaces, on the other hand, there seems to 
be no satisfactory notion of "locality" . To circumvent this problem, we focus on 
observables instead. An observable A E 13(H) is localized in a region A if A acts 
like the identity on all tensor factors outside of A. 

Reversible dynamics on such a system is again represented by a unitary 
U € U(H). We use the symbol a to denote the action by conjugation of U 
on observables: A n- a(A) = UAU* . 

Let {cr x j Cy^ j Cz } be the Pauli matrices acting on the ith spin. Since one can 
form a basis of B(H) from products of the Pauli matrices {cr x ,cr z }%L\ acting 
on single spins alone, it suffices to specify the effect of a on these 2N matrices 
in order to completely determine the dynamics. For example, we can set 

a(o-^)=at 1) ®^ ) ®4 l+1 \ (4) 
=*•«. (5) 

It is a simple exercise to verify that the operators on the right hand side fulfill 
the same commutation relations as the {o4 , <Xz }. This is sufficient to ensure 
that a unitary U implementing the time evolution a actually exists. 

As in the previous example, we can ask more refined questions about U. 
For example: can we implement U by a fixed-depth circuit of nearest-neighbor 
unitaries? Can one interpolate between U and the trivial time evolution while 
keeping observables localized along the path? Is there a simple numerical invari- 
ant which would allow us to easily decide these questions? 

In this particular case an educated guess gives rise to affirmative answers to 
all these questions. Indeed, set 

/l 
(M+1) _ 1 
Vt "0010 
\0 e i4 

with respect to the standard basis of the ith and (i + l)th spin. Define 

N 

Ut =J{v^ +1) (7) 

i=l 

(the product is well-defined because Vt commutes with its translates, so the order 
in which the product is taken is immaterial). Then one checks that a corresponds 
to Ui, whereas Uq is the trivial evolution. Clearly XJ t continuously interpolates 
between these two cases, and, by construction, does not increase the localization 
region of an observable by more than two sites. We note that U\ is known in the 
quantum information literature as the interaction used to generate graph states 

H5B21- 

In Section [3 wc will set up a general theory for answering the questions 
posed above — including in cases where one is not so lucky to have an explicit 
parametrization in terms of gates at hand. 
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3. Systems 

3.1. Quantum Walks. We consider a quantum system with a spatial degree of 
freedom igZ, and at every site or "cell" x a finite dimensional "one-cell Hilbcrt 
space" H x ■ These spaces can be thought of as describing the internal states of 
the system as opposed to the external, spatial variables x. The Hilbcrt space of 
the system is 

oo 

H= 0%, (8) 

X — — OQ 

The specification of the 'Hx will be called a cell structure. 

We call a unitary operator U on % causal^, or a quantum walk, if, for any 
i€Z, there are some X- < x + such that <j> £ T-L x C H implies Ucj) £ ®y= X - 
We assume that x± are both non-decreasing as functions of x, and go to ±oo 
when x does. 

By Uy X ■ Hx Hy w e denote the block matrix corresponding to the direct 
sum (|SJ|, i.e. 

^®^ = ©EM- (9) 

x y x 

The causality of U implies that, for any x, only finitely many y give non-zero 
summands. 

Grouping. — The spatial variable a; of a walk and the internal degrees of freedom, 
described in H x , arc largely interchangeable. In one direction wc can choose a 
basis \x, 1), . . . , \x, al) in some H x , and replace the point x by the sequence of 
points (a;, 1), ... , (x, d), each with a one-dimensional space H. x ,i of internal states. 
Because in the above definition, we assumed that the spatial variable ranges over 
Z, groupings have to be followed up by a relabeling of sites in the obvious way. 
In the other direction, we can "fuse together" several cells xi, . . . , Xk, getting a 
new cell X with internal state space Hx = Hi- I n either case it is clear 

how to adjust the neighborhood parameters x±. 

Hence we can either regard our system as one without internal degrees of 
freedom, and Hilbert space % = £ 2 (1). Typically this may involve some large 
neighborhoods [x—, x+\. Or else, we can group cells until we get a nearest neigh- 
borhood system, i.e., x± = x ± 1, at the expense of having to deal with high- 
dimensional H x - 

Most of these definitions and constructions are easily generalized to higher- 
dimensional lattices. It is therefore instructive to identify the feature which re- 
stricts our results to the one-dimensional case. Indeed, it lies in the fact that 
one can choose a partitioning into intervals [a, , bi] C Z such that the sites below 
<ij and above 6j interact only through the interval [a,,6i]. Formally, that is a 
consequence of demanding lim a; _i.± 00 x± = ±oo. It is easy to see that such a 

We are indebted to a referee who drew our attention to the inflationary use of the word 
"local" in our manuscript, where, among other things, quantum walks were called "local" 
unitaries. We changed this to "causal" for the crucial finite propagation property of walks and 
cellular automata. This is reminiscent of relativistic propagation in algebraic quantum field 
theory and in keeping with usage in signal processing (where x would be time, and X— = x). 
In quantum information it agrees with [4J2]. The terminology disagrees with what some field 
theorists would probably say 1181 . and with 1161 . 
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separation need not be possible in a two-dimensional lattice, even if every cell 
has a finite neighborhood. In this more general setup, neighborhood relations up 
to regroupings may be described in terms of coarse geometry [TJ], a theme we 
will not pursue here. 



Translationally invariant walks. — The simplest way to define a cell structure 
is to choose a Hilbert space "Ho, and to set T~L X = Hq for all igZ. We then have 
the unitary equivalence T~L = ^ 2 (Z) <g) Hq. In that case we can define the shift 
operation S and its powers by 

S n (\x}®cf>o) = \x + n)®<f> . (10) 

In this setting one frequently looks at translationally invariant walks, i.e., uni- 
taries U commuting with S. More generally, there might be some period p such 
that [U, S p ] = 0. Clearly, it is natural in this case to group p consecutive cells, 
so that after grouping one gets a strictly translationally invariant walk. 

The space T-Lq can then either be considered as an internal degree of freedom 
of a walking particle, or as a coin so that shift steps (possibly depending on 
the internal state) are alternated with unitary coin tosses I <g> C. Translationally 
invariant systems will be treated in more detail in Sect. 16.41 



Periodic boundary conditions. — Since we are after a local theory of quantum 
walks, global aspects — like the distinction between a walk on Z and a walk 
on a large ring of M sites — are secondary, as long as the interaction length 
L = max|a;_|_ — X—\ remains small in comparison with M. In fact, from any 
walk on a ring we can construct one on Z which locally looks the same. More 
formally, let the sites of the ring be labeled by the classes Zm of integers modulo 
M, and identified with {0, . . . , M — 1}. Then we extend the cell structure by 
setting Hx+kM = H x for all k € Z. In order to extend the unitary U on the ring 
to a walk U on Z we set 

ft _ / if \x-y\ > L , , 

xy ~ \ U x , y , if \x-y\< L, x' = x,y'=y modA/. 

For the second line to be unambiguous, we require 2L < M. To verify unitarity 
we need to compute 

J2uA = Y, u *'y' u y'*'> ( 12 ) 

y y' 

where we have assumed that \x — z\ < 2L, because otherwise the left hand side 
is zero anyhow. Note that for each summation index y only one class y' € Zm 
can occur in the sum on the right hand side. Moreover, every class y' appears, 
although possibly with a zero contribution. But the sum on the right is S X ' Z '^x', 
which together with a similar argument for UU* proves the unitarity of U. 

From the point of view of index theory, the walk U on the ring and U on the 
line are the same. However, if we iterate U, the interaction length (x + — X—) 
increases, and eventually non-zero matrix elements can occur anywhere in U n . 
In this sense, the set of quantum walks on a ring does not form a group. This 
is the reason why the theory of walks on Z is more elegant and more complete. 
From now on we will therefore consider walks on Z only. 
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3.2. Cellular Automata. Once again we consider a system in which a finite di- 
mensional Hilbert space H x is associated with every site ieZ. However, rather 
than combining these in a direct sum, we take their tensor product. In plain 
English this means that, for any two sites x,y, rather than having a system of 
type H x at position x or a system of type T-L y at position y, as in a quantum 
walk, we now have a system of type H x at position x and a system of type Hy 
at position y. In contrast to the infinite direct sum of Hilbert spaces, the infinite 
tensor product is not well-defined. Since we want to look at local properties, we 
could work with a "potentially" infinite product, i.e., some finite product with 
more factors added as needed in the course of an argument. But it is easier to 
work instead with the observable algebras A x , equal to the operators on the 
Hilbert space H x , or equivalcntly the algebra M.d( x )! where d(x) = dimH x {Aid 
denotes the algebra of d x d-matrices) . In analogy to the definition for walks, we 
will refer to the specification of the algebras A x as the cell structure. 

For the observable algebras associated to sets A C Z, we use the following 
notations: for finite A, A{A) is the tensor product of all A x with x € A. For 
Ai C A 2 we identify A(A X ) with the subalgebra A(A{) ® l A2 \ Al c A(A 2 ). For 
infinite A C Z we denote by A{A) the C*-closure of the increasing family of finite 
dimensional algebras A(Af) for finite Af C A, also called the quasi- local algebra 
[5]. In particular, the algebra of the whole chain is «4(Z), sometimes abbreviated 
to A. 

A cellular automaton with cell structure {A x } x ^z is an automorphism a of 
A = A(Z)) such that, for some functions x H> x± as specified in Sect. 13.11 
each a(A x ) C *4([x_ , x + ]). Note that the restricted homomorphisms a x : A x —> 
A([x-,x+]) uniquely determine a, because every observable acting on a finite 
number of cells can be expanded into products of one-site observables. These local 
rules a x have to satisfy the constraint that the algebras a x {A x ) for different x 
commute element-wise. In that case they uniquely determine an endomorphism 
a. For examples and various construction methods for cellular automata wc refer 
to [16]. 

Exactly as in the case of quantum walks we can group cells together for 
convenience. Whereas the dimensions for subcells add up for quantum walks 
(dim(£) xeA ri x = J2 x eA dimHx) we get A(A) = Md with the product d = 

By considering the time evolution of observables, we have implicitly chosen to 
work in the Heisenbcrg picture. The expectation value of the physical procedure 
(i) prepare a state p, (ii) run the automaton for k time steps, (m) measure 
an observable A would thus be given by the expression p(a k (A)). Accordingly, 
we choose a convention for the shift on a chain with isomorphic cells, which 
at first seems inverted relative to the definition ([TU)) for walks. We define it 
as the automorphism a with a(A x ) = A x -\, acting according to the assumed 
isomorphism of all the cell algebras. Thus if one prepares a certain state, it will 
be found shifted to the right after one step of a, in accordance with ([TU]) although 
in that case UT-L X = T-L x +\. 

4. Local Implementability 

We have defined the causality properties of walks and cellular automata ax- 
iomatically, i.e., as a condition on the input-output behavior of the maps U and 
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a. Alternatively, one may take a constructive approach. Here, one would list a 
set of operations that should certainly be included in the set of local dynamics, 
and refer to any given time evolution as being locally implementable if it can be 
represented as sequence of these basic building blocks. Both methods are equally 
valid, and in this section we will completely analyze their relation. From the ax- 
iomatic point of view this might be called a "structure theorem" , whereas from 
the constructive point of view one would call it a "characterization theorem" . 

In the case at hand, there is a natural choice of building blocks. Namely, 
we can partition the system into some subsets ("blocks") of sites, and apply 
a unitary operation separately to each subsystem in the partition. (Note that 
the unitaries would be combined by a direct sum for walks and by a tensor 
product for cellular automata). For such maps the interplay between unitarity 
and causality is trivial: causality puts no constraint whatsoever on the choice 
of unitaries acting in each block. Moreover, it allows the overall operation to 
be resolved into a sequence of steps, in which one block operation is done after 
the other. This picture is close to the gate model of quantum computation |13j : 
here each block unitary would correspond to one "gate" involving some subset 
of registers, so that these gates do not disturb each other. The fact that they can 
be executed in parallel is expressed by saying that these infinitely many gates 
nevertheless represent an operation of logical depth 1. 

For partitioned unitary operations the various block unitaries obviously com- 
mute. Commutation is really the essential feature if we want to resolve the overall 
time step into a sequence of block unitary steps. Indeed, consider a family of 
commuting unitaries Uj , each localized in a finite subset Aj of some infinite lat- 
tice (not necessarily one-dimensional). We only need that the cover by the Aj is 
locally finite, i.e., each point x is contained in at most finitely many Aj. Then 
the product FT . Uj implements a well-defined operation on localized elements. 
In the cellular automaton case (where localization just means Uj € A(Aj)), we 
define the action on a local observable A as 

a(A) = (j[U*)A(j[U j ) > (13) 

j j 

with the understanding that both products range over the same index set, namely 
those j for which Aj meets the localization region of A. Here the products can 
be taken without regard to operator ordering, since we assumed that the Uj 
commute. Including additional factors Uj on both sides will not change a(A), 
since such factors can be "commuted past" all other Uj> , Up and A to meet 
the corresponding U* and cancel. So the product is over all j, in the sense of 
a product over any sufficiently large finite set. Similar considerations apply for 
the case of walks on general lattices. 

Now if a QCA a is represented in the form (fT3"|) , we can also represent it as 
a product of partitioned operations: indeed, we only need to group the Uj into 
families within which all Aj are disjoint. The product of each family is obviously 
a partitioned unitary and under suitable uniformity conditions on the cover we 
only need a finite product of such partitioned operations to represent a, typically 
s + 1 factors, where s is the spatial dimension of the lattice. Hence we consider 
the representation as a product of partitioned unitaries as essentially equivalent 
to the representation by commuting unitaries as in (|13[) . In cither case we will 
say that the system is locally implementable. 
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We now come to the basic result for implementing general walks or cellular 
automata by commuting unitaries — provided we are allowed to enlarge the 
system. The key feature of these extensions is that they work in arbitrary (not 
necessarily one-dimensional) lattices and that the ancillary system is a copy of 
the system itself, on which we implement simultaneously the inverse operation. In 
the following result, we allow the underlying "lattice" X to be any countable set. 
General neighborhood schemes are described as controlled sets in some coarse 
structure . For the present paper it suffices to describe causality in terms of a 
metric d on X, of which we only assume that all balls Ml{x) = {y\d(x, y) < L] 
are finite sets. The causality condition for walks on X is then that there is some 
"interaction radius" L such that in (|9|) U yx = for d(x,y) > L. Similarly, for 
QCAs, the causality condition reads a(A x ) C A(Nl{x)). For walks U,V on 
the same lattice we simply write U © V for a walk with one-cell Hilbert spaces 
7-tx®1C x , where W x are the one-cell spaces for U and JC X those for V. This splits 
the total Hilbert space into "H®JC, and U ® V acts according to this direct sum. 
Similarly, we define the tensor product a © /3 acting on two parallel systems 
combined in a tensor product. 

Proposition 1 

(1) For any quantum walk U , the walk U © U* is locally implementable. 

(2) For any cellular automaton a, the automaton a © a -1 is locally imple- 
mentable. 

Proof. (1) We are considering a doubled system in which the one-cell Hilbert 
space at x is 'H. X ®'U. X . Let S x denote the unitary operator on the doubled system 
which swaps these two summands, and acts as the identity on the one-cell spaces 
of all other sites. Now consider the unitaries 

T X = (U*®1)S X (U®1). 

These commute, because they are the images of the commuting transformations 
S x under the same unitary conjugation. Moreover, they are localized near x by 
the causality properties postulated for U. Hence their infinite product defines a 
walk unitary, as discussed above. This unitary is 

~[ T x = (U* © 1)5(17 © 1) = 5(1 @U*)(U@ 1) = S{U © IT), 

X 

where we have used that S = Y[ x S x is just the global swap of the two system 
copies. Hence U © U* = (Y\ x S X )(Y\ X T x ) is locally implemented. 
(2) Essentially the same idea works for cellular automata Again we consider 
the unitaries S x £ A x ®A X , which swap the two tensor factors, so that S X (A X © 
B x ) = (B x © A X )S X . Now consider the unitary elements 

T x = (id® a) [5*]. 

Here we have written the arguments of an automorphism in brackets, to dis- 
tinguish it from grouping parentheses, and thus eliminate a possible source of 
confusion in the coming computations. As images of a family of commuting uni- 
taries under an automorphism, the T x arc themselves a commuting family of 
unitaries. Moreover, they are localized in A x © Aj^t x y Hence they implement 
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a cellular automaton /?. We determine it by letting it act first on a localized 
element of the form A x £g> 1 with A x <G A x . 

p[A x ® 1] = (J[ T„) * (A x ® 1) (H Ty) 

v y 

= (id g> a) [[J SJ (id ® a) [a* (8) l] (id ® a) [jj S* y ] 

= (id <g> a)[l g> A x ] = 1 ® a[A x ] 

A similar computation shows that (3[l ® = -B x (g) 1. Since a is an au- 

tomorphism, this is the same as /3[1 ® Bx] = or 1 !- 8 *] ® I. Using the ho- 
momorphism property of j3, we get for general localized elements A, B that 
f3[A<& B] = oT l [B] <g) a [A]. Hence following f3 by a global swap (implemented by 
[\ x S x ) we have implemented a® a -1 locally. □ 

5. The group of locally computable invariants 

In this section we take up the idea of a locally computable invariant and show 
that, for cither walks or automata, these invariants necessarily form an abelian 
group. The group multiplication reflects both the composition and the parallel 
application to a double chain. We postpone to later sections the question whether 
nontrivial invariants exist, i.e., at this stage it might well be that the group 
described here is trivial. Later on we will determine this group to be (Z, +) 
for quantum walks (see Sect. |5])and (Q+, •), the multiplicative group of positive 
fractions, for cellular automata (see Sect. [7]). In this section, in order not to 
double each step, we will describe the arguments for the case of walks, and just 
comment at the end on the necessary changes for the cellular automaton case. 

Suppose we have defined a property P(U), which is defined for any quantum 
walk U, and which can be determined solely on the basis of a finite collection 
of the block matrices U xy . More specifically, if we write the walk in nearest- 
neighbor form by grouping, we call a property P(U) locally computable if we can 
compute it from the restriction of U on any interval of length > 2. The crucial 
part of this definition is, of course, that the result obtained in this way must be 
the same for any interval we may select for the computation, a property which 
we stress by calling P a locally computable invariant. 

Suppose now that two walks U\ and [7 2 share a patch, in the sense that there 
is a long interval [xi,X2] C Z, on which the Hilbcrt spaces 1-L X for x £ [x\,X2\ 
have the same dimensions and, after the choice of a suitable isomorphism, the 
unitaries U\ and U2 restricted to these subspaces act in the same way. We assume 
that the interval is sufficiently long to determine P. Then local computability just 
means that we must have ~P{U\) = P(U-2)- In other words, P must be constant 
on each equivalence class of the relation of "sharing a patch" . So the theory of 
locally computable invariants is really equivalent to characterizing the classes of 
the transitive hull of this relation: We will write U\ ~ U n , if there is a chain of 
walks Ux, U2, ■ ■ ■ , U n such that, for all i, Ui and f7j+i share a patch. In contrast 
to the relation of sharing a patch, this equivalence relation no longer makes 
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Fig. 2. Combining a crossover from U\ to Ui with a crossover from Ui to t/3 to obtain a 
crossover from U\ to U3. The shaded double chains can be fused to a single cell by Prop, IT1 



any requirements about the sizes of any one-cell Hilbert spaces in the walks Ui 
and U n . The most comprehensive locally computable property is now just the 
property of U to belong to some equivalence class: all other locally computable 
properties are functions of this class. Our aim thus shifts to computing the set J 
of equivalence classes for "~" . The equivalence class of a walk U will be denoted 
by ind (U) € J, and called its (abstract) index. 



Triviality of locally implementable systems. — Let us first make the connection 
to the questions of the previous section: suppose that a walk or automaton 
is locally implementable, i.e., the product of a collection of block partitioned 
unitaries. Compare this with a system in which all unitaries, whose localization 
intersects the positive half axis, are replaced by the identity. Clearly, this acts 
like the identity on all cells on the positive axis, and we can further modify the 
system by making it trivial (O-dimensional T~i x or 1-dimensional algebras A x ) 
for x > 0. Clearly these systems share a large patch (most of the negative axis), 
so they are equivalent. In other words, locally implementable systems have the 
same index as the identity on a trivial chain. 



Crossovers. — A very useful fact about the relation U\ ~ U2 is that it implies a 
prima facie much stronger relation: It is equivalent to the property that there is 
a "crossover" walk U c , which coincides with U\ on a negative half line {x\x < a} 
and coincides with U2 on a positive half line {x\x > b}. 

Proof. Let us denote the relation just described by U\ rj U-2- Then Ui rj U% => 
U\ ~ JJ2, because U\ and U2 each share an infinite patch with U c . 

In the converse direction, if U\ and U2 share a patch, we can define U c to 
be the walk whose one-cell Hilbert spaces H x are those of U\ for x to the left 
of the shared patch, and are those of U2 for x to the right of the shared patch. 
Similarly, we define the unitary U c to coincide with U\ to the left and with U2 
on the right. On the shared patch we can choose either one, since the two walks 
coincide. Since the shared patch was assumed to be sufficiently long this does 
not lead to an ambiguity for either U c or U~ x . Hence U\ w U2. 

In order to cover the case that Ui and U2 arc linked by a chain in which 
any neighbors share a patch, we only need to prove that U± ~ U2 is a transitive 
relation. 
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In order to prove transitivity, consider the walks Ui,U2,U3 with crossovers 
U12 and U23 as shown in Fig. [2 There we also included a copy of the inverse of 
\Ji- We take the overall picture as a representation of the walk C/ 12 © U^ 1 © t^23- 
Now consider the strands of U12 and U^ 1 to the right of the crossover region 
of U\2- Since U2 © U^ 1 is locally implementable by Prop. [TJ we can replace this 
pair of strands by a trivial system, still retaining a legitimate unitary operation 
for the rest. Similarly, we can fuse the strands of U23 and U^ 1 to the left of 
the crossover region of U23 to nothing. This results in a unitary operator, which 
coincides with U\ on a left half axis and with U3 on a positive half axis, i.e., a 
crossover U13. □ 

Grouping. — Suppose we regroup some finite collection of the cells. Clearly, 
this does not affect cells far away, and we immediately get a crossover. Hence 
the index does not change when regrouping cells, even if this is carried out in 
parallel. We will implicitly use this in the sequel by regrouping sites in whatever 
way is most convenient. 

Parallel chains. — Define the direct sum U © V of walks as in the previous 
section. Then if U c (resp. V c ) is a crossover between U\ and U% (resp. V\ and 
V2), U C @V C is obviously a crossover between U\ © Vi and U2 © Vi- Hence the 
class ind (U © V) € J depends only on the equivalence classes of U and V, and 
we can define an "addition" of indices by ind (U) + ind (V) = ind (U © V). This 
addition is abelian, because there is a trivial crossover between U © VandV © U, 
just exchanging the summands on a half chain. Moreover, since an inverse is 
defined as — ind(J7) = ind([/*) via Prop. Q] we conclude that J becomes an 
abelian group. 

Products. — Now suppose U and V are walks on the same cell structure, so that 
UV makes sense. We claim that there is a crossover between UV and U © V. 
Hence we also get 

ind (UV) = ind (U © V) = ind (U) + ind (V). 

Indeed, consider the cell structure on which U © V is defined, which has one- 
cell Hilbert spaces Hx © V. x at each site. Now let S + denote the unitary which 
acts as the swap on all T-L x © T-L x with x > and leaves the subspaces with 
x < unchanged. Then (U © 1)5*+ (1 © V)S+ is a crossover between U © V and 
UV®1~UV. 

Modifications for cellular automata. — The concept of crossovers and the argu- 
ments for the group structure can be taken over verbatim, with the replacements 
U 1 — y OL) Hx !— ^ A x , © H- ®. Of course, the index group will also be different, and 
we will adopt the convention to write it multiplicatively. The product formula 
thus reads ind (a/3) = ind (a <£>/?) = ind (a)ind (f3). 

Numerical representation and shift subgroups. — To make the abstract theory 
of this section useful, one needs to establish an isomorphism of the index group 
with some explicitly known group. The natural way to do that is to identify 
generators, i.e., some particular walks which cannot be implemented locally, but 
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are sufficient to generate arbitrary walks up to locally implcmcntable factors. 
Although it is far from obvious at this point, it will turn out later that, for walks 
as well as for cellular automata, the role of generators is played by the shifts. 
Since there is only one kind of generators, it suffices to assign numbers as "index 
values" to the shifts to establish an isomorphism of the abstract index groups 
with groups of numbers. 

For walks, the shift was introduced in Sect. 13.11 We denote by Sd the shift 
on a system with c?-dimensional internal degree of freedom (dim H x = d for all 
x). Similarly, let o~d denote the shift automorphism on a chain with cell algebra 
Aid- We tentatively demand 

ind {Sd) = d (14) 

and, similarly, 

ind (a d ) = d. (15) 

This has to be consistent for shifts on parallel chains. Since Sd © S e = Sd+ e and 
<Jd®a e = &d-e, this requires that we take the indices of walks as a group of num- 
bers under addition, and for the cellular automata as a group of numbers under 
multiplication. Indeed, we will show that the above formulas fix an isomorphism 
of the abstract index group J to (Z, +) for quantum walks, and to the group 
(Q+, •) of positive fractions for cellular automata. 



6. Index for quantum walks 

6.1. Pedestrian definition. The following is the basic definition of this chapter. 
To the best knowledge of the authors, it is due to Kitaev [12], who calls this 
quantity the flow of a walk U . 

Definition 2 For any walk U , we define the index as 

indC/= (^(U xy yU xy -TT{U yx )*U yx ). (16) 

x>0>y 

Note that the sum is finite by virtue of the definition of causal unitaries. Clearly, 
for the simple shift we get ind Si = 1, confirming Eq. (TT4"|) . 

Of course, we will show presently that this quantity has all the properties of 
the abstract index discussed in the previous sections. However, from the defini- 
tion given here it seems miraculous that such a quantity should be always an 
integer, and independent of the positioning of the cut. To see this it is better to 
rewrite this quantity in the following way. 



6.2. Operator theoretic definition. We introduce the projection P for the half 
axis {x > 0}, i.e., the projection onto the subspace @ x>0 T~L X . Then, for 4> x £ % x 
and 4> y 6 H y we get 

x > and y > 

-{<t>x\U xy \(j) y ) x < and y > , 
((t>x\U xy \4> v ) x>0andy<0 
x < and y < 



\PU-UP\4> y ) = 
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Hence, for any pair (x, y) the commutator [U, P] has just the signs used in the 
definition of the index, and we get 



Note that for the the trace on the right hand side we cannot use linearity of the 
trace to write it as the difference of two (equal!) terms, because this would result 
in an indeterminate expression oo — oo. 



6.3. Fundamental properties of the index for walks. 
Theorem 3 

1. vaAU is an integer for any walk U 

2. ind U is locally computable, and uniquely parameterizes the equivalence classes 
for the relation ~ from Sect. [5l hence can be identified with the abstract index 
defined there. 

3. ind {U\ © U2) = ind (Z7i) + ind (U2), and, when U\ and U2 are defined on the 
same cell structure, ind (U1U2) = ind (U±) + ind (U?). Moreover, for the shift 
of d-dimensional cells: ind Sd = d. 

4. ind U = if and only if U admits a "local decoupling", i.e. there is a uni- 
tary V, which acts like the identity on all but finitely many H x , such that 
UV is block diagonal with respect to the decomposition % = (@ x<0 Hx) © 



case it can be written as a product of just two partitioned unitaries. When U 
is regrouped in nearest neighbor form, then the partitioned unitaries can be 
chosen to couple only pairs of nearest neighbors. 
6. mAUi — ind Uq if and only ifUo and U\ lie in the same connected component, 
i.e., there is a norm continuous path [0,1] 3 t 1— > Ut of causal unitaries of 
uniformly bounded interaction length L with the specified boundary values. 

The rest of this subsection is devoted to the proof of this result. According 
to Eq. psp. the index is closely related to a difference of projections. If these 
were finite dimensional, we could just use linearity to get the difference of two 
integers. The following Lemma shows that the result is still an integer when the 
difference of the two projections has finite rank. Actually it is even sufficient for 
the difference to be trace class, and with a careful discussion of the trace, it is 
sufficient for the ±1 eigenspaces of the difference to be finite dimensional [3]. 
Here we include the simple case sufficient for our purposes. 

Lemma 4 Let Q, P be orthogonal projections in a Hilbert space 1-1, such that 
Q — P has finite rank. Then 

1. The range 1Z = (Q — P)H is an invariant subspace for both Q and P. 

2. Tr(Q — P) is an integer. 

3. There is a unitary operator V such that V(f> = </> for all (f> _L (Q — P)H, and 



such that Q > V*PV orQ< V*PV. 
4- //Tr(Q — P) = 0, the V from the previous item satisfies Q = V*PV. 



ind U = Tr U* [P, U] = Tr (U*PU - P) 



(18) 



5. 
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Proof. 1. follows from the identity 

Q(Q -P) = Q(i -P) = (Q- P)(i - P), 

and its analogue for P. 

2. Clearly, we can evaluate the trace in a basis of TZ since the basis elements from 
TZ 1 - contribute only zeros. Since the restrictions of Q and P to TZ are projections 
on a finite dimensional space, 

Tr(Q - P) = Tr K (Q - P) = Tr K (Q) - Tr K (P) 

is the difference of two natural numbers. 

3&4. Obviously, we can find such a unitary on TZ with the corresponding 
property for the restrictions of Q and P to TZ. We then extend V to be the 
identity on TZ . When Tr(Q — P) = 0, this V is a unitary mapping from QTZ to 
PTZ. a 

Proof (of Theorem^). None of the statements, or values of the index will change 
under grouping, except part of item 5, which requires nearest neighbor form. 
Therefore we will assume without loss that all walks are nearest neighbor. We 
will use Lemma |H with Q = U*PU. 

1. This follows directly from Eq. ([T5]) and Lemma HI item 2. 

2. Let P' be the projection onto another half axis, say x > x$. Then P — P' 
is finite rank and hence (U*PU - P) - (U*P'U - P') = P' - P - U*(P - P')U 
is the difference of two finite rank operators with equal trace. Hence the index 
docs not depend on the cut position, and since formula ^ clearly involves only 
matrix elements at most 1 site from the cut, it is a locally computable invariant. 
It remains to be shown that it is a complete invariant, i.e., that ind U\ = ind U 2 
implies U\ ~ U 2 in the sense of Sect. [5] This will be done in connection with 
item 4 below. 

3. This follows from Sect. [5] But a direct proof (for the product) is also 
instructive: 

(u 1 u 2 )*P(u 1 u 2 ) -p = (u;pu 2 -p) + u 2 *{u* 1 pu 1 - p)u 2 

is the sum of two finite rank operators, of which we can take the trace separately. 

4. Apply Lemma QJ item 3, to get V with P = V*(U*PU)V, and hence 
PUV = UVP. The fact that V — 1 vanishes on all but finitely many Ti x follows 
from its construction: V—l vanishes on the complement of (P—Q)Ti C W-iSiTio, 
for the cut "-1|0" used in Dcf.H Note that this implies V ~ 1, and also UV - I, 
since a unitary which has no matrix elements connecting x > and x < clearly 
allows a crossover with the identity. From the product formula for locally com- 
putable invariants we therefore get that ind U — implies U ~ I. Obviously, this 
extends to other values of the index: if ind[/i = ind U 2 we have ind U*U 2 = 0, 
hence UIU 2 ~ 1 and hence U\ ~ U 2 . This completes the proof of item 2. 

5&6. These items each contain a trivial direction: We have already shown 
in Sect. [5] that locally implcmcntablc walks have trivial index. Moreover, it is 
clear from Definition [2] that the index is a continuous function, and must hence 
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be constant on each connected component. The non-trivial statement in 5. is 
that walks with trivial index are indeed implementable, and in 6. that walks 
with vanishing index can be connected to the identity (the rest then follows by 
multiplication). In either case, an explicit construction is required, and it will 
actually be the same one. 

So let U = U\ be a walk ind U = 0. Let Vq denote the decoupling unitary for 
the cut — 1|0, obtained in the proof of item 4, and define similar unitaries Vk for 
the cuts at 2fc— l|2fc such that UVk has no non-zero matrix elements (UVk) xy with 
y < 2k < x. Let Hk denote a hcrmitian operator located on the same subspaces 
as Vk — I, such that Vk = exp(iHk). We will take all Hk bounded in norm by the 
same constant (ir will do). Then since they live on orthogonal subspaces, their 
sum H = X^fc-^fc i s well-defined and also bounded. Now let V(t) = cxp(itH), 
which is a norm continuous function of t, because \\H\\ < ir. The endpoint V(l) 
can also be defined by this product formula V(l) = Ilfc^' because on each 
subspace %2k-\ ® %2fc only one these unitaries is different from 1. Moreover, 
UV(1) has no matrix elements y < 2k < x for any k, i.e., it is block diagonal for a 
decomposition of Z into pairs {2k, 2k + 1}. Now take a similar Hamiltonian path 
deforming each block in this matrix decomposition separately to the identity. 
Specifically, we take W(0) = 1 and W{1) = UV{1). Then t i-> W(t)V(t)* is a 
norm continuous path (although no longer a unitary group), with the endpoints 
1 and U. Moreover, each unitary W(t) or V(t)* is based on a partition into 
neighboring pairs so that, for no t, W(t)V(t)* has any non-zero matrix element 
between sites with \x — y\ > 2. This proves the remaining statement in item 5 
(for t = 1), and also the statement about uniformly bounded neighborhoods in 
item 6. □ 



6.4- The translation invariant case. Suppose that U commutes with some power 
of the shift. It is then useful to group spaces T-L x into larger blocks to get com- 
mutation with the shift itself. That is, in this section we assume all H x = JC to 
be equal, and U xy = U x ~ y , where by a slight abuse of notation the single-index 
quantity U x is defined as U x q. The width L is the largest x such that U x ^ 
or U- x ^ 0. It is natural to diagonalizc U using the Fourier transform. We de- 
fine T : £ 2 {Z) ® JC -> £ 2 ([-7r,7r]) ® K, by F{^){p) = e lpx V(x). This is 

to be read as a /C-valued equation, where we use the natural identification of 
£ 2 ([— 7r, 7r]) £g> JC with the set of /C-valued square intcgrable functions on [—it, n]. 
Similarly, we identify £ 2 (Z) (g) K. with the /C-valued square summable sequences. 
Then J-U T* becomes the multiplication operator by the p-dependent matrix 

L 

Uip) = £ U x e**. (19) 

X — — L 

Note that this is a Laurent polynomial in e lp . The largest degree of e lp in the 
polynomial is x + — x, which is constant by translation invariance. The lowest 
degree is x — X- . Further, U(j>) must be a unitary operator on K, for every 
p € M. Taking these facts together, we conclude that both the determinant 
detU(p) = f(p) and its inverse l/f(p) = det U(p)* are Laurent polynomials as 
well. But this is only possible if / is actually a monomial, say proportional to 
exp(mp) for some integer n. We claim that this n is the index: 
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Proposition 5 For a translation invariant walk, 

dctU{p) = Ce lpind( - u \ (20) 

for some phase constant C. 

As a simple example consider the shift on a chain with dim "Ho = 1- We already 
noted after Definition [2] that this has index 1. The corresponding p-dependent 
unitary is the number U(p) = e %v , so this also gives index 1. For unitaries acting 
on each site separately in the same way, we get agreement because ind U = 0, and 
U (p) is independent of p. Note also that both sides of [2D] have the same behavior 
under composition and direct sums. This proves the formula for all walks which 
can be composed of shifts and sitcwisc rotations. Actually, all translationally 
invariant walks can be represented in this way [7] , but we prefer to give a direct 
proof of the proposition without invoking this decomposition. 

Proof. From Definition [2] we get 

oo —1 

ind U = J2 J2 \U x - y \ 2 -Tr \U y - x \ 2 

x—0 y—~oo 
oo 

= J2 nTv\U n \ 2 =Y,S nm nTiiU^U n ) 

71— — co nm 

On the other hand, for any invertible matrix function U, 

4-detU(p)=detU(p) TYfuip)- 1 ^^-). 
dp V dp 1 

Hence with det U(j>) = exp(ipn) the above integrand is actually constant equal 
to in, and ind (U) = n. □ 

The properties of a walk depend crucially on the properties of the eigenvalues 
e^W, . . . , e lUd ^ of U(p) (see Fig.©. Clearly, p i-» U(p) is an analytic family of 
operators, so we can follow the eigenvalues as analytic curves [IT]. The derivatives 
of the eigenvalues determine the group velocity: Let P a (p) be the eigenprojection 
of U (p) belonging to the eigenvalue exp iuj a (p) . Then the group velocity operator 
can be written as 

G( P ) = lira ijT(t) = £p q(p )%M, (21) 

t-s-oo t *■ — ' dp 
a 

where X denotes the position operator, and at degenerate eigenvalues the projec- 
tions P a (p) have to be chosen appropriately, as dictated by analytic perturbation 
theory. The limit is in the strong resolvent sense. Hence the probability distribu- 
tion for the selfadjoint operator G in a state p is equal to the asymptotic position 
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Fig. 3. Example for eigenvalues of U(p), with dim/C = 3, L = 5, ind!7 = —1. In this case the 
spectrum is a single curve on the torus. The index can be computed by the signed number of 
crossings of any horizontal line, or as the sum of the derivatives of all branches. 

distribution starting from p in "ballistic" scalingjl]. In particular, when the in- 
ternal state is unpolarised, i.e., when the initial state is of the form p = a® I/d, 
we find P2 

(X(t)) = {X(0))+ t -ind(U). 

In this sense the index has direct relevance as a mean speed of the walk. 

It is, of course, suggestive to connect the topological properties of the eigen- 
value curves in Fig. [3] with the index. In principle, these curves are allowed to 
cross each other. So in general, we have several branches of curves, which wind 
several times around the torus before closing. The number of curves and their 
winding numbers would appear to be a topological invariant, but this is not true 
because of the "avoided crossing" phenomenon, in which a small perturbation 
of the walk turns an eigenvalue crossing into a close encounter of two separate 
curves (suggested also at p = 1.6 in Fig. [3]). Hence only the sum of the winding 
numbers is stable with respect to small perturbations, and this is indeed equal to 
ind U. To sec this, note that at every p the sum of the derivatives of all branches 
is equal to the index. Therefore the sum of the winding angles of the branches 
taken from p — — 7r to p = n is 27rind U. The winding angles of the closed curve 
components are just suitable sub-sums of this, and this partition is changed at 
avoided crossing points. 

We close this section by establishing a variant of the main theorem for the 
translationally invariant case. Of course, most of this follows by simply specializ- 
ing. The only question which requires a new argument is whether the path con- 
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necting two translationally invariant walks can be chosen to go entirely through 
translationally invariant walks. This is established in the following proposition. 

Proposition 6 Let U be a translationally invariant causal unitary with ind U = 
0. Then we can find a norm continuous path tt-}Ut of translationally invariant 
causal unitaries of bounded width such that Uq = 1 and U\ = U . 

Proof. We use the factorization theorem for "paraunitary" operators [7] (see also 
[T7]V which states that U(p) can be written as a finite product 



U(p) = V ]JW mk (p)V k 



k=l 



of constant unitaries and the elementary causal unitaries 



W m (p) = 



' 1 



• ox 

• ' 



V ■■• 1/ 



As a quantum walk these corresponds to a shift of only the first internal state by 
m positions, leaving all other internal states fixed. This walk has index m, and 
the product formula yields indC/ = X)fc TO fc- The maximal polynomial degree of 
any matrix element in e ±tp depends on the localization region [x+, —X-], and is 
clearly bounded by £ max = ^2 k \rrik\- Hence we can contract the walk to 1 by 
contracting each of the Vk to I, never exceeding width L max on the way. □ 



7. The index for cellular automata 

For cellular automata we proceed in analogy to the case of walks, i.e., by defining 
directly a locally computable quantity as the index ind a of a walk automorphism 
a. We then establish that it is actually a complete locally computable invariant 
and, at the same time that it characterizes the connected components of QCAs. 

As a technical preparation we need some background on commutation proper- 
ties of algebras spanning several cells. It is basically taken from [IB], and included 
here to make the presentation here self-contained. 



7.1. Support algebras. For defining the index we need to find a quantitative 
characterization of "how much" of one cell ends up in another. To this end we 
introduce the notion of support algebras. Consider a subalgcbra A C B\ <g>2?2 of a 
tensor product. What is the position of A relative to the tensor structure? Here 
we answer a relatively simple part of this question: which elements of B\ , B2 are 
actually needed to build Al 

For the following Definition with Lemma, recall that A' denotes the commu- 
tant {a|Vai S A : [a, a{\ = 0}. 

Lemma 7 Let B± and B2 be finite dimensional C*-algebras, and A C B\ ® Bi 
a subalgebra. Then 
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1. There is a smallest C*-subalgebra C\ C B\ such that A C C± ® Bi- We call 
this the support algebra of A on B\, and denote it by C\ = S(A,B\). 

2. Consider a basis {e^} C B2, so that every a € A has a unique expansion 
a = ^2 a M ® e M with e 61. Then S(A, B\ ) is generated by all the elements 
dp arising in this way. 

3. The commutant of S(A,Bi) in B\ is characterized as {b € B\ \ b® 1 € A'}. 

Proof. We can pick out the terms a M by applying a suitable functional u; M from 
the dual basis to the second factor, i.e., by applying the map id ® w M : B\ ® B2 — > 
B\, which takes b\ ® 62 to ^(62)61- Clearly, if a £ C\ ® B2, so that a can be 
expanded into simple tensors using only elements from C\ in the first factor, we 
find = (id® w M )(a) G C\. Hence the algebra described in item 2 must be 
contained in any C\ satisfying item 1. Since it also satisfies condition 1, wc have 
identified the unique smallest C\. The characterization 3 follows by looking at 
commutators of the form [b ® 1, a], and expanding a as above. □ 

This construction was introduced to the QI community by Zanardi |19j , who 
applied it to the algebra generated by an interaction Hamiltonian, and conse- 
quently called it an "interaction algebra" . Of course, we can apply the construc- 
tion also to the second factor, so that 

AcS(A,B 1 )®S(A,B 2 )cB 1 ®B 2 . (22) 

The crucial fact we need about support algebras is that "commutation of alge- 
bras with overlapping localization happens on the intersection" . More precisely, 
we have the following 

Lemma 8 Let Ai C B\ ®B2 and A2 C B2 ® £>3 be subalgebras such that A\ ® I3 
and li ® A2 commute in B\ ® B2 ® B3 . TTien S(-Aij $2) andS(A2,B2) commute 
in B2- 

Proof. Pick bases {e M } C B\ and {e' v } C B2, and let a G Ai and a' € Ai. 
Then we may expand uniquely: a = e M ® a M and a' = a! v ® e' v . Then by 
assumption 

= [a <g> 1 3 , Ii ® a 1 } = ® K> of v ] ® e' v . 

Now since the elements e M ® e' v are a basis of B\ ® B3, this expansion is unique, 
so we must have [a^,a^] = for all /i, v. Clearly, this property also transfers to 
the algebras generated by the a M and a' v , i.e., to the support algebras introduced 
in the Lemma. □ 



7.2. Defining the Index. Now consider a cellular automaton a on a cell structure 
A x - By regrouping, if necessary, we may assume that it has only nearest neighbor 
interactions. Now consider any two neighboring cells A2X ® A2X+1, and their 
image under a, i.e., 

a(A 2x ® ^2x+i) C (Azx-i ® A 2x ) ® (A2 X+ i (8 A 2x +2) • 
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Fig. 4. Cell structure with support algebras. (Read top to bottom) If the width of cells is 
taken as log dimension, the index can be read off the slant in the boxes representing mapping 
by a. 



We apply the support algebra construction to this inclusion, obtaining two al- 
gebras 



K 2x = s(a(A 2x ®A 2x+ i), {A 2x -i ® A 2x )) 
K 2x+ i = S ( a(A 2x ® A 2x+ i) , (A 2x+ i <g> A 2x + 2 ] 



(23) 
(24) 



Intuitively, the algebras lZ y with even index become larger when information 
flows to the right, whereas the ones with odd index describe a flow to the left. 
Of course, this intuition will be made precise below, but at this stage one can 
at least check these statements for the shift: When a(A y ) = A y -\ (recall the 
convention made at the end of Sect. I3.2[) we have lZ 2x = A 2x -i <8> A 2x and 



TZ- 2x = CI. This is to be contrasted with lZ y = A y for the identity. 

Continuing with our construction, observe that by Lemma [SJ the subalgebras 
1Z 2x +i and lZ 2x + 2 commute in the algebra (A 2x +i ® A 2x + 2 ) . Algebras TZ y which 
are further away commute anyhow, since they arc contained in disjoint cells. We 
conclude that all lZ y commute. 

By definition of support algebras, a(A 2x <8> A 2x +i) C lZ 2x (g> TZ 2x +i, so that 
the algebras TZ X together generate an algebra containing aA(Z). Because a is 
an automorphism, this is the same as _4(Z). Now if any 1Z X had a non-trivial 
center, i.e., if there were an element X £ 1Z X commuting with all of 1Z X but not a 
multiple of 1, this X would also be in the center of the entire quasi-local algebra 
A(Z). However, this center is known to be trivial [5]. We conclude that each 
1Z X must have trivial center, and hence be isomorphic to M. r ^ x ) for some integer 
r(x). Moreover, the inclusion noted at the beginning of this paragraph cannot 
be strict, since otherwise we would find an element in the relative commutant, 
which would once again be in the center aA(Z). To summarize, we must have 



hence 



t(A 2x <g> A 2x ^ 
d{2x)d{2x -f 



— Tl 2x <8> Ti 2x +i, 
= r(2x)r(2x + 1) 



(25) 
(26) 



On the other hand, the commuting full matrix algebras T^x+i and lZ 2x + 2 to- 
gether span the tensor product isomorphic to M r { 2x +i)r( 2x + 2 ) inside (A 2x +i <E> 
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•A2X+2) ■ Again the inclusion cannot be strict, because otherwise the automor- 
phism would not be onto. From this we get the second relation and dimension 
equation 

K2X+1 ® K 2x +2 = A 2x +i ® A 2x+2, (27) 
hence r(2x + l)r(2x + 2) = d{2x + l)d(2x + 2). 

These relations arc summarized pictorially in Fig. |4l They give us the the first 
two equalities in 

r(2x) d(2x + l) r(2x + 2) _ . 

d(2x) = H2x-TT) = df^TT) = md a - (28) 

In other words, this quantity is constant along the chain, and hence we are 
entitled to define it as the index ind a. The even/odd asymmetry only comes from 
the construction, by which the even lZ 2x describe a flow to the right (increasing 
x), and the odd IZ2X+1 are associated with a flow to the left. By shifting the 
entire construction, we could switch the even/odd distinction, and define, for 
any A x , both the ascending and the descending 1Z X . In any case, for the shift ad 
of a d-dimensional chain we get d(y) = d, r(2x) = d 2 , and r(2x — 1) = 1, and 
hence ind a c i — d, as announced in Eq. (jTSj) . 

It is part of the local computability property that we have a lot of freedom 
in choosing the cell structure for which we want to evaluate the index. Since one 
typically wants to use this freedom to simplify the computation, we will now 
summarize the constraints. It is clear that there are three subalgebras involved 
in the computation, playing the role of, say the above Aq, Ai, A 2 for determining 
TZ\. Let us call these Al , Am , Ar to emphasize that these algebras need not be 
part of the original cell structure, and we are free to choose them within certain 
limits. Let us start by fixing some algebra Am — Md, a full matrix algebra 
contained in some local algebra, whose crucial property is to split the system: 
we must have an isomorphism of the total algebra with A<m ® Am <8> A>m, 
where the outer factors contain suitable infinite half chain algebras, and such 
that 

[a(A <M ),A >M ] = [a(A >M ),A <M ] = {0}. 

Clearly, this imposes a a lower bound on the size of Am m terms of the inter- 
action length of the automorphism. Now we choose finite dimensional matrix 
subalgebras Ar C A>m and Al C A<m such that ci(Am ) C Al ® Am 65 Ar. 
These three can be taken as part of a nearest neighbor cell structure, so that 
the above arguments give 

Am = M d \ . , d 



S(a(AL®A M ),A M ®A R ) = M, 



=>inda = - (29) 



r 



Note that there is no harm in choosing Ar larger than necessary: the support 
algebra, being the minimal algebra needed to build the tensor product, will 
simply not change. By a similar argument, we can choose Al too large without 
changing this support algebra. 
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7.3. Fundamental properties of the index for cellular automata. 
Theorem 9 

1. inda is a positive rational for every a. When the automaton is regrouped in 
nearest neighbor form, both the numerator and the denominator of ind a in 
canceled form divide every cell dimension. 

2. ind a is locally computable, and uniquely characterizes the equivalence classes 
for the relation ~ from Sect. [3J It can hence be identified with the abstract 
index defined there. 

3. ind (a ® a') = ind (a)ind (a 1 ) and, when a and a' are defined on the same 
cell structure, ind(aa') = ind (a) ind (a'). 

If, for some y, a(A((— oo, y])) C A((— oo, y])), then inda £ N. 
Moreover, for the shift of d- dimensional cells: hidcr^ = d. 

4- An automaton a has index 1 if and only if it can be implemented locally. 
In this case it can be written as a product of just two partitioned unitary 
automorphisms. If a is partitioned in nearest neighbor form, the partitioned 
automorphisms can be taken to couple pairs of nearest neighbors only. 

5. Two automata 0,0, 0,1 on the same cell structure have the same index if and 
only if they can be deformed to each other, i.e., there is a strongly continuous 
path [0, 1] 3 t t— > at of automorphisms, all with the same neighborhoods, and 
with the specified boundary values. 

Proof. (1,3) inda € Q+ follows immediately from the construction in Sect. 17. 21 
particularly Eq. (|28|) . Let inda = | be the fraction in canceled form. Then, from 

this equation p = nr(2x) and q — nd(2x), where n is the canceled factor. Hence 
q divides d(2x) and from the second fraction in (|28p we find that p divides d(2x + 
1). By shifting the construction by one cell, we find the remaining divisibility 
statements. From Eq. ([25]) we also get the product formulas 3. Suppose that a 
maps some left half chain in to itself. Then we choose a cell partition so that 
y = 2x + 1 in the setting of ([23")) . Then H2X+1 C A2X+1 and since these are full 
matrix algebras with the same unit, the quotient inda = dix+\ / l T2x+\ in (|28|) is 
integer. The value for the shift was verified as an example after (|28[) . 

(2,4,5) From Eq. ([25)1 local computability is obvious. From the general discus- 
sion in Sect.[5]we also get that locally implementable a have inda = 1. The hard 
part, which is needed to identify the abstract index with the concrete formula is 
the converse. For this the crucial step is the following claim: 

Let a and a' be nearest neighbor cellular automata on the same cell structure 
and with the same index. Then there are unitaries V x € A x <£> A x +i such that 
the two locally implemented automorphisms 

f3{A) = (Y[V 2x . 1 rA([[V 2x . 1 ) 

X X 

7^)=(n^x)*^(n^x) 

X X 

satisfy: a'7 = /3a. 

Before proving this claim, let us see how it implies the statements in the 
Theorem. By the general theory of Sect. [5] locally implementable automorphisms 
/3, 7 are ~-equivalent to the identity, and a ~ /3a = a'7 ~ a'. Hence equality 
of the indices as defined by (|2T))) implies crossover equivalence, and hence the 
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equality of all locally computable invariants. This proves item 2. The converse 
in item 4 follows by taking a 1 = id, giving the local implementation a = /3~ lr y 
of any index 1 automorphism a. 

Finally, it is clear for item 5 that we can connect a and a' with the same 
index by the required continuous path: we just need to contract each unitary 
V x in /3,7 to the identity, to obtain a path at = fita^ 1 with ao = a and 
ot\ = a'. This path will not be continuous in the norm on automorphisms, 
i.e., we cannot make \\a t — a s \\ small, since this would already fail for one-site 

operations a t (A) = (uf°°^j AU®°° with t h> [/ ( norm continuous. However, 

for any finitely localized observable A, t t— > at (A) will be continuous in norm, 
which is the claim of strong continuity made in the Theorem. We remark that 
an important part of the proof of item 5 is missing at this point: We did not 
exclude the possibility that there are continuous paths linking automorphisms 
of different index. This will be achieved by Prop. 1131 an expression for the index 
which is manifestly continuous with respect to strongly continuous deformations. 

Now to prove the claim, let a and a' have the same cell structure and the same 
index. Then we carry out the construction of Sect. [7721 for both automorphisms, 
resulting in some intermediate algebras TZ X and TV X . 

Since the indices coincide, formula (|28[) demands that these are full matrix 
algebras of the same dimensions. For example, 

TZ 2x -i ® TZ 2x = A 2x -i ® A 2x = R-2X-1 ® TV 2x . 

Clearly, there is a unitary operator V 2x -\ £ A 2x -\®A 2x so that V 2x _{R y V 2x -i = 
TV y for y = 2x and for y = 2x — 1. We can take all these unitaries together as 
implementing one partitioned automorphism (3(A) = (Y[ x V2x-i)*A([[ x V 2x -i). 
By definition, it satisfies (3(H y ) = TV for all i/gZ. 

Now consider the action of a and a' on A 2x (8 A 2x +i- We now get two iso- 
morphisms 

a' : A 2x ® A 2x +\ -t TV 2x <g> TV 2x+l and 
(3a : A 2x ® A 2x+ i -t TZ 2x ® TZ 2x+1 -> TV 2x ® TV 2x + x . 

Hence /3a restricts to an automorphism of A 2x ®.A2a;+ij and can therefore 

be implemented by a unitary V 2x £ A 2x <8> A 2x +i- These unitaries together 
implement 7, and we get the desired equation a'j = (3a. □ 



7.4- Index for classical reversible automata. In this section, we will review a 
common notion of index for reversible classical cellular automata, and show 
that it coincides with our definition. In the context of this paper, a reversible 
classical cellular automaton can be defined as a particular case of a quantum 
cellular automaton. In each cell A x we single out a maximal abelian subalgebra 
V x . With respect to a suitable choice of basis, T> x is then the set of diagonal 
matrices. As a finite dimensional abelian algebra, we can regard T> x as the set 
of complex valued functions on a finite set A^, called the alphabet of the cell, 
which at the same time serves as the set of basis labels for the orthonormal basis 
in which T> x is diagonal. The global C*-algebra of the classical system is then 
the infinite tensor product P(Z) = (^^L.^ Z^c- It is canonically isomorphic 
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to the algebra of continuous functions on the compact cartesian product space 
Az = y^L_ oc A x , also known as the space of infinite configurations. We use 
lower case letters such as c for such configurations, and denote by c(x) G A^ the 
configuration of the cell at x. 

Now let a be a QCA with the property that a(T>) C T>. Then the restriction of 
a to I? is an automorphism of T>, which must be of the form (a/)(c) = /(^(c)), 
where <P is a homcomorphism on configurations. The causality conditions on a 
are readily expressed in terms of and, together with analogous arguments 
for the inverses show that $ is a reversible classical cellular automaton in the 
usual sense, apart from the requirement of translation invariance. There are some 
subtle points to note about the correspondence a — > <P: 

— Suppose that in the above argument we start from a general, not necessarily 
strictly causal automorphism a of the quasi- local algebra A(Z). We still get a 
continuous <P on the compact space Ag. By the definition of the product topol- 
ogy this means that the local configurations <&{c) x after the time step depend 
on only finitely many c{y). In the translation invariant case this means that 
is a cellular automaton with finite neighborhood. In fact, this argument is 
used to establish that the inverse of a reversibly cellular automaton also has 
a finite neighborhood. This is a rather surprising sharpening of the causality 
condition. However, we are appealing here to a highly non-constructive com- 
pactness argument, which gives no control on the size of the neighborhoods, 
or (barring translation invariance) on the uniformity of the neighborhoods. 
For example, we can apply it to a cellular automaton a on a 2D lattice, whose 
quasi-local C*-algebra is isomorphic to that of a ID automaton. Hence the 
condition of a being an automorphism is not strong enough to give a ID 
automaton in the sense defined above. 

— The mapping is onto, i.e., every classical reversible cellular automaton can 
be "quantized". The argument is very simple for finite lattices, e.g., a regular 
lattice with periodic boundary conditions: one labels the basis of a Hilbcrt 
space by the classical configurations. Then the classical automaton <P is a 
permutation of the basis vectors, which can be interpreted as a unitary oper- 
ator via U$\c) = |^(c)). Then for all observables we set a(A) = U^AU^,. One 
needs to check that this transformation is causal in the quantum sense |16j . 
in particular that off-diagonal local operators that are finitely localized (i.e. 
localized on a finite number of cells) keep this property under the action of 
a. Indeed one gets a bound on the quantum neighborhoods, which involves 
both the neighborhoods of <P and the neighborhoods of <P~ 1 . The same com- 
putation provides a formula for a(A), for A finitely localized, in terms of the 
classical rule, and this can be used to define a also for infinite lattices. We 
will not, however, make this explicit here. 

— The mapping a — > <1> is not injective. Indeed an ambiguity is inherent in 
the construction just described: we can choose different bases with the same 
diagonal operators \a)(a\ by choosing a phase for each basis vector. This 
amounts to changing a by an on-site unitary, which is certainly irrelevant for 
index purposes. But we can consider this more generally: suppose that two 
cellular automata a, /3 restrict to the same automorphism on the diagonal 
algebra T>. Then a/3 -1 leaves V point-wise fixed. In a finite lattice, so that a 
is unitarily implemented, the implementing unitary hence commutes with all 
elements of the maximally abelian algebra T>, hence must itself be diagonal. 
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It is suggestive that this also holds in a localized form on the infinite lattice, 
so a/3^ 1 would be a product of commuting unitaries and hence have trivial 
index. The following result builds on this intuition. 

Proposition 10 Suppose that a and [3 be quantum cellular automata taking the 
diagonal algebra T> into itself, with the same restriction to T>. Then ind a = ind (3. 

Before coming to the proof of this proposition we single out two arguments 
of independent interest, each of which can be used to draw the main conclusion, 
without discussing in detail the structure of local phase factors. The first criterion 
uses the absence of propagation. The second uses the global transpose map 
: A{1j) — >• Z. It is defined as the matrix transpose on each local algebra, in a 
basis in which T> is diagonal. Since transposition is isometric, and consistent with 
the embcddings A i— >• A <8> 1 it extends to the whole algebra. is a linear anti- 
homomorphism (moaning 0(AB) = 0(B)0(A)), and, for every automorphism 
a, 0a0 is again an automorphism. 

Lemma 11 1. Let a be a nearest neighbor cellular automaton such that, for 
some finite interval [z_, z + ] we havea n (A([0, 1])) C A([z-,z+]) for all n £ Z. 
Then ind a = 1 . 

2. For any cellular automaton a: ind (0a0) = ind a. 

Proof. 1 . The index ind a™ can be expressed as a ratio of of subcell dimensions 
of „4([z_, z+]), hence is uniformly bounded in n. But since ind a™ = (ind a)™ this 
implies ind a < 1. With the same argument for the inverse we get ind a > 1. 

2. By assumption, the global transposition is made with respect to product 
basis, so that for a tensor product B = B\ ® #2 of cells we get 0q = 0® x ® 
0B 2 - We will drop the indices on in the sequel. Then it is clear from the 
definition (Lemma [7]) that the support algebra construction behaves naturally 
under global transposition, i.e., we have S(0A, 0B\) = 0S(A, Bi). Moreover, 
when B\ is a tensor product of cells we get 0B\ = B\ . Hence in (f2"5| we find for 
the automorphism a — 0a0: 



Since this has the same dimension as "T^-, we find from (|28p that ind (0a0) = 
ind a. □ 

Proof (of Prop. \TU\) . Due to the multiplication formula, we only need to consider 
the case that a is equal to the identity (= (3) on T>. Consider any operator 
A x <G A x in a single cell, and let D y £ T> y be a diagonal element in another cell. 
Then 

a{A x )D y = a(A x )a(D y ) = a{A x D y ) = a(D y A x D y ) = D y a(A x ). 

Hence the finitely localized element a(A x ) commutes with all diagonal operators 
on the neighboring sites, and a(A x ) £ A x ®(£) , x T> y . This algebra is best seen as 
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a direct sum of copies of A x , labeled by configurations c = {c(y)} with c(y) £ A y 
of all cells y ^= x in the localization region of a(A x ). A homomorphism of A x 
into this algebra splits into one homomorphism into each summand, which in 
turn is given by a unitary U x (c). Hence we can summarize the action of a on 



a{A x ) = U x (c)*A x U x (c) <g> P(c), (30) 

c 

where P c is the minimal projection of the diagonal algebra corresponding to c. 
We also know that diagonal elements of A x are fixed, so U x (c) is itself diagonal, 
say U x (c)\a) — u(a, c)\a). This leads to 

a((|o>(6|)x) =EtM (\°)(b\)*®P(c)- (31) 

The commutation of a(A x ) and a^-Ay) for x ^ y introduces further conditions 
on the phase functions u. But rather than analyzing these in detail, we use the 
Lemma to conclude directly from Eq. (f3~T|) that inda = 1. To this end, note that 
by applying the homomorphism a to (|31|) and using a(P(c)) = -P(c), we get a 
corresponding formula for the iterate of a: 

«"((l«)(&l)x)-E(^|) n (\a){b\) x ®P(c), (32) 

for any n £ Z. Clearly, the localization region of this operator does not increase 
with n, so by the first part of the Lemma we get inda = 1. Alternatively, we can 
apply to the equation, using 0(P(c) = P(c). This reverses each of the ketbra 
operators, so OaO = a~ 1 . Hence inda = 1 also follows with the second part of 
the Lemma. □ 

An index for classical reversible cellular automata has been defined albeit only 
in the translationally invariant case. According to G. A. Hedlund ([5], section 14), 
the definition is due to L. R. Welch, so we will call it the Welch index i\y here. For 
the definition itself we will follow Kari [TO], where it is introduced in section 3. 
We will show that this coincides with the quantum index. Hence the quantum 
index is a possible extension to non-translationally invariant systems. It is very 
likely that the theory in [10] can also be extended directly, but we have not gone 
to the trouble of checking all the details. 

In the translation invariant case, all cell alphabets A x = A are the same, and 
a cellular automaton is a map <f> : A z — > A z . Let r be a "large enough" integer, 
and the set of 4r-tuples of the form 

(c(0), • • • , c(2r - 1), (0c)(-r), • • • , (<t>c)(r - 1)) , 

where c runs over all infinite configurations. Then the Welch index of <fi is 
defined as 

iw(4>) = S . (33) 
|A| 

Clearly, for the identity only |A| 3 different letters occur here, so i^(id) = 1. 
Similarly, for a shift we get i\y(S) = |A|, and parallel application of cj> and 4>' to 
parallel chains yields iw{4> x 4>') = 'wf^'ifW)- 
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The non-trivial results about the index, and the structure theory of reversible, 
translationally invariant classical cellular automata is developed in |10) , with key 
results analogous to our paper: the expression ([55)1 does not depend on r (pro- 
vided it is large enough). The product formula holds, and an automaton <fi is 
locally implementable iff iw{4 > ) = !• Moreover, every such automaton can be 
decomposed into shifts and locally implementable ones. Since a classical "local 
implementation" implies a partitioned representation of the quantum automa- 
ton, we can put these facts together to conclude that 



for any quantum cellular automaton which restricts on the diagonal subalgebra 
to a classical CA given by (f>. In this sense our theory is a direct generalization 
of Kari's work, extended by the aspects of deformation classes (which make no 
sense in the classical discrete setting) and local computability (which makes no 
sense in the translation invariant setting) . 

7.5. Interlude: More Analogies between Walks and Cellular Automata. The two 
definitions of the index, ([T5|) for walks and (|2l)]) for cellular automata are not 
directly analogous. Here we would like to point out the differences, and discuss 
how to make the analogy between these two cases even tighter by supplying the 
missing analogous definitions. 

The definition (f!2T)| considers a part of the system split into three parts L-M-R. 
Based on suitable inclusions, it gives a formula for the index, which immediately 
makes obvious that it is always a positive rational. In contrast, the walk expres- 
sion Eq. (|16[) is a difference of numbers which can take arbitrary real positive 
values, and only one cut of the system is considered. Moreover, (|16p made it very 
easy to prove the continuity of the index under deformations, whereas neither 
(l29l) nor the abstract considerations of Sects. 2] and [5] clarify continuity for the 
index of cellular automata. Since continuity is an important feature of our index 
theory, we will need an appropriate expression also for the automaton case, and 
the analogies laid out in this subsection are intended to motivate the form of 
this formula. 

Dimension based formula for walks. — Let us first set up an index formula for 
walks in analogy with (|29p . The analog of the support algebra is the "support 
subspace of K-12 in £2", denoted S^K-^, K-2), which is defined for any subspace 
K,\2 C JC\ © IC2 of the orthogonal direct sum of Hilbert spaces. Namely, it is the 
smallest subspace C C K.% such that K.12 C K\ © C Then the analog of Lemma[5] 
holds in the sense that subspaces K.12 C K\ © K2 and K23 C IC2 © K3 are 
orthogonal iff S(/Ci2,/C2) -L ^(^23,^2) arc orthogonal. Now consider subspaces 
T4l © T~Lm © Hr C H chosen with the localization constraints as in Sect. 17.21 In 
particular, we require UHl -L Hr, and UHr _L Hl- Then the direct analog of 
(ggjl reads 



Here the second equality, in which Plm is the projection onto %l © Hm etc., 
follows with S(U(H L ®n M ),n M © Ur) = PlmUPmbH- 



iwii 1 ) = hid a, 



(34) 



indU = &voiS{U(U l ®Um),11m®'Hr) -&\mU M 
= rank (PmrUPlm) - dim% M - 



(35) 
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Proof (of Eq. \35\) ). Consider the block matrix for U with respect to the decom- 
position 

where the pieces at the ends contain the appropriate infinite half chains. Using 
the causality assumptions, we find that 



U 



/ * * 





* * * 





* * 


* 


0* 




loo 


* * 



(36) 



where the asterisks stand for any possibly non-zero block. We have highlighted in 
boldface the block PmrUPlm appearing in (|3"5|) . and introduced two separating 
lines, namely the separation — oo, L, M\R, +oo on the domain side (i.e., a vertical 
line) and the horizontal separation — oo, L\M, R, +oo on the range side of U. 
These separators do not cross on the diagonal, which is why we cannot simply 
compute the index from the highlighted block via Eq. (|16p . However, this is 
easily amended by multiplying with a suitable shift: we can introduce a basis 
in each T-L x , and hence in H, effectively making all underlying cell dimensions 
one-dimensional. In this representation we can introduce a shift operation S, and 
clearly U' = S^ M ^U will be a unitary with the same matrix elements as U shifted 
vertically by \M \ = dim'Hju. Obviously, hid U' = hid f/+|M|, which explains the 
second term in (|35[) . and leaves us with proving that ind U' = rank {PmrUPlm)- 
Clearly, for this task the further subdivision of the blocks is irrelevant, and we 
can consider a general block decomposed unitary operator 



U 








U 2 i 


U22 



with a hnitc rank upper right corner. We have to show that rank U21 = Ti^L^i^i)- 
But from the unitarity equation it follows that U21 is a partial isometry, so 
U21U21 is its domain projection, whose dimension is indeed the rank of Ui\- d 



Half neighborhoods and one-cut dimension formulas. — The proof of formula 
(pi5j) was essentially by reduction to the case of "half neighborhoods", i.e., the 
case that = [x,x + 1], in which no influence ever spreads to the left. 

Then one of the off-diagonal blocks of the unitary entering (|T6| vanishes, and 
we saw directly that the other block gives an integer contribution, which can be 
interpreted as a dimension. 

Similarly, for a half- neighborhood cellular automaton we demand a(A x ) C 
Ax— 1 <8 A x , where once again we have chosen the convention to match informa- 
tion traveling to the right, observing the Heisenberg picture. In particular, this 
condition is satisfied by the shift. For such automata with we can simplify the 
index formula in a way quite analogous to the case of half-neighborhood walks. 
Indeed, setting 

% = S{a(A x ),A x ) nndM x = S(a(A x ),A x -x), (37) 

we can employ the same arguments as in Sect.[7]to conclude that these commute, 
and must be full matrix algebras T x = Mt(x) an d Af x — ■M n (x)- Then, further 
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following the previous reasoning, a(A x ) = Af x <E> T x and A x = T x Af x +i- This 
yields the dimension equation 

d{x) = n(x)t(x) = t(x)n(x - 1). (38) 

Hence the integers n(x) do not depend on x. In fact, hid a = n(x). This follows 
readily from the observation that 

K 2x +i = s(a(A 2x <£> A 2x +i), (A 2x +i ® A 2x+2 ) ) = % x+ i O J2Z+2- 

Therefore, from Eq. we get 

inda = d(2x + l)/r{2x + 1) = d(2x + l)/t(2a; + 1) = n(2x+ 1). 

In fact, for any cell structure on which a shift is available, we could have 
used this slightly simplified setup to define the index of any a by first shifting 
and regrouping to a half-neighborhood automaton, and correcting by a factor 
depending on the size of the necessary shift. However, since such a shift is not 
available in general, there was no gain in this approach. 

For nearest neighbor automata the structure of support algebras using just a 
single cut, i.e., algebras of the form S(a(^4i), Ar) are not sufficient to define the 
index. As a counterexample, consider a unitary evolution A t— > U* AU for U G 
Ah ® Ar with Ah = Ar = Aid- Clearly, as a locally implementable operation, 
this has always trivial index. But the support algebra written above can be just 
about anything. For example, taking a "controlled unitary" U = J2i ® Ui: 
we have S((x(Ar),Al) as the algebra of diagonal matrices, and S(a(Ah), Ar) 
generated by the U*Uj, which for Ar = Ai r can easily be Ar. On the other 
hand, one can easily construct an automorphism with the same support algebras, 
but index d. 

The algebraic structure of support algebras across a single cut is also insuf- 
ficient for us in another way: it is a discrete structure, hence does not go to 
any trivial value as we deform an a to the identity. Therefore, if only to get the 
continuity of the index we need to look at some continuously varying quantities. 
We found a formula with just these properties by looking at cellular automaton 
analogs of the simple walk formula (fTo| . The following section is devoted to the 
proof of this formula. 



7.6. One-cut quotient formula. Notation for normalized traces. — A continu- 
ously varying quantity measuring the difference between subalgebras and going 
to a trivial value as they come to coincide may be some kind of angle, or overlap, 
between the linear subspaces. It is natural to measure such angles with respect 
to the only scalar product between algebra elements, which is canonically de- 
fined in our context. Indeed, let t denote the normalized trace on the algebra 
A of the entire chain. That is, on any matrix subalgebra B C A, B = Aid, we 
define t(A) = ^ Tr(^4), where Tr is the usual matrix trace on Aid, which is 1 on 
minimal projections. The reason for this normalization is that in contrast to the 
matrix trace, the value of r does not change if we consider A to be embedded 
as A Cg) 1 in some larger subalgebra B <g> B±, and is hence a well defined state on 
chain algebra A. In fact, it is the unique state (normalized positive functional) 
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on A with the property that t(AB) — t(BA). This characterization is purely 
algebraic, which implies that, for every automorphism a of A and any A £ A, 
we have r{a{A)) = t(A). 

We now define the scalar product between algebra elements by 

(x\y) = r(x*y). (39) 

The completion of the algebra A as a Hilbcrt space with this scalar product 
is called the GNS- Hilbcrt space H T of the tracial state r. We write \x) € W T 
for the vector obtained by embedding x € A in the completion. The trace of 
operators on T-L T will be denoted by Tr, to avoid confusion with the matrix trace 
Tr of some elements of A, which is also used later. Since r is preserved by any 
automorphism a, we can define a unitary operator V a on H T with 

V a \x) = \a(x)). (40) 

Consider now a finite or infinite dimensional subalgcbra B C A. By P we 
denote the orthogonal projection onto the closure of B in H T . For finite matrix 
algebras B = M.d the matrix units = clearly form a basis B, and one 

readily verifies that {Vd \eij)}f - =1 is an orthonormal basis of this d 2 -dimensional 
subspace. It is sometimes also helpful to represent the projection P as an integral 
over unitaries. That is, for a finite dimensional subalgebra B C A we have 

P = d 2 JdU \U)(U\, (41) 

where the integral is over the unitary group of B, and dU denotes the normal- 
ized Haar measure. For infinite dimensional subalgebras such formulas are not 
available. We will only need infinite dimensional projections of this type for half 
chain algebras, which are the closure of an increasing net of finite dimensional 
subalgebras B n . In this case the family of projections P n associated with the 
approximating algebras is also increasing, and converges strongly to P. 

When B\ C A and B2 C A are commuting matrix subalgebras, the corre- 
sponding matrix units, say ej^ and e^ b , can be taken together as a set of matrix 
units for B\®Bi, and we get 

(e {1) \e {2) ) -rfe (1) e (2) ) - —6 -6 , 

Therefore we get 

PiP 2 = |i)(l|. (42) 

This equation also holds for infinite pieces of the chain, e.g., a right and a left 
half chain localized on disjoint subsets of Z. This readily seen by approximat- 
ing each half chain by finite matrix algebras and using the strong convergence 
of projections. Note that (|42p also implies that Pi and P2 commute, and the 
projections Pi — |J)(J| with the intersection removed are orthogonal. 

Overlap of algebras. — Of course, if two algebras do not commute, which requires 
that the localization regions of Bi and B2 overlap, the geometric position of the 
subspaces PiH T is not so simple. Even if the algebras have trivial intersection 
£11162 = CI, so that (PiH T ) n {PiH-r) = C|I) the vectors in the remainder 
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can now have angles different from tt/2, and may even approximate each other. 
This leads to the following definition of a quantitative notion of the overlap of 
two algebras, which will be the basis of the index formula we develop in this 
section. We state it together with a few of its basic properties. A matrix algebra 
here always means a subalgebra, which is isomorphic to Md for some d < oo 
and contains the identity of A. 

Lemma 12 For any to subalgebras B\,B2 C A, with corresponding orthogonal 
projections P±,P 2 on T-L T , we define their overlap as 



V(B U B 2 ) = ^Tr(P 1 P 2 ) g [0,oo]. 



(43) 



Then 



1. r)(B\,B 2 ) = 1, for commuting matrix algebras 

2. 77 (Si , i?2 ) > 1> for any two subalgebras (with unit). 

3. Let B\,B 2 C B12 and B%,Bi C $34 are matrix subalgebras such that B12 and 
B34 commute. Then 



n(B 1 B 3 ,B 2 B i ) = r)(B 1 ,B 2 )r)(B 3 ,B4 



(44) 



Let 0i, B2, B^, B4 be a chain of matrix algebras such that [Bi, Bf\ = {0}, except 
for i,j = 2,3. Then v(BiB 2 ,B 3 B 4 ) = f](B 2 ,B 3 ). 

Let £>i, B 2 , B3 be commuting matrix algebras with B2 = Md- Then 77(61 £>2, B2B3) 
d. 

Let B\ = Md cind B2 be finite dimensional matrix algebras, and a an auto- 
morphism of the ambient algebra such that ||a(a;) — x\\ < e||a;|| for x € B\. 
Then 

\n(a(B 1 ),B2)~v(Bi,B2)\<ed 2 . 



6. 



Proof. Item 1 follows directly from Eq. (|42]) . For item 2, note that Pi > |J)(1|. 
Since we are only considering subalgebras containing the identity of the ambient 
algebra, the parenthetical remark is only added for emphasis. Item 3 follows 
by observing that the normalized trace on B12B34 = B\ 2 <£> S34 is fixed to be 
the product of normalized traces on the subalgebras. In this tensor product 
representation one readily verifies that the projection belonging to B1B3 is the 
tensor product Pi and similarly for P 2 ,Pi. Then the formula follows because 
the trace TV also factorizes. 

Finally item 4 follows by taking B1B4 and B2B3 as the pairing of item 2, and 
using item 1 to conclude that T}(fB\,B?) = 1. 

To prove the continuity estimate in item 6, note that overlaps are > 1, so 
that by a gradient estimate on the square root function 



\r){a(B 1 ),B 2 )-r){B 1 ,B 2 )\ < - Tr (( P ± - P 1 )P 2 ) 



- 2" 



Pi\ 



where P\ 1 P\,P 2 are the projections belonging to a(Bi) 7 Z?i, B2, and ||^||i de- 
notes the trace norm of X. We use the representation of Pi in the form (|41[) . so 
that 



Pi - Pi 



IdU 



(\a(U))(a(U)\-\U)(U\ 
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Wc will estimate the trace norm of this expression by estimating the integrand, 
and using that, for any unit vectors <p, ip in a Hilbert space, || \<j>) (<j>\ — \tp) < 
2||</> — Indeed, for unitaries like U and a{U) the vectors \U) and \a{U)) 

have norm 1 in U T . Moreover, \\\a(U) - U)\\ 2 = r((a(U) - U)*{a(U) - ufj < 

e 2 ||[/||r(l) = e 2 . Hence 

l\\Pi-Pi\\i<^e JdU = d 2 e. 

□ 

Index formula. — With these preparations we can state the main result of this 
section: 

Proposition 13 Let A< = -A(-oo.o] and Ay =-A[i,oo)> and Al-, Ar the algebras 
of any two neighboring cell algebras for a nearest neighbor grouping of the chain. 
Then 

7?(a(.A<),.A>) r)(a{A L ),A R ) 

ind (a) = -) ( = -) (. (45) 

i 1 [a(A > ),A < ) T][a{A R ),A L j 

Moreover, if ' t i— ^ a t is a strongly continuous family of cellular automata with the 
same cell structure and neighborhood scheme, then ind (at) is constant. 

Proof. The square of the numerator of the second index expression is TtV^PlV^Pr, 
and we first verify that this expression is unchanged if we choose a larger L and R, 
say L' = LiUL and R' = RURi. Indeed, the algebras a{Ah x ), u(Al), Ar, Ar x 
satisfy the conditions of Lemma 1121 item 4. Moreover, arguing as for formula 
(l42"T) we see that not only the trace, but even the operator is independent of an 
enlargement. Hence taking a strong limit we obtain the corresponding expression 
for the infinite half-chains. 

From Lemma [T2l item 3, it is clear that the expression thus defined satisfies 
the tensor product property in fact, numerator and denominator do so indepen- 
dently Moreover, the formula is valid for shift automorphisms by virtue of item 
5 of the Lemma. 

Now let a be a tensor product of shifts with index = ind a, and let ind 'a be 
the value the formula in the proposition gives for any automorphism a. Then 
we have ind'(a ® c" 1 ) = ind '(a) (ind cr) _1 = ind 'a/ (ind a). So it remains to 
prove that ind 'a = 1 for every a with ind a = 1, which by Thm. [9] means that 
a is implemented by two layers of block unitaries. Moreover, we can forget all 
unitaries acting only on one side of the separation, since they do not change the 
algebras. Only one unitary Ulr € Al <£> Ar connecting the cells immediately 
adjacent to the cut remains. Hence it only remains to prove the Lemma below. 

For the continuity statement observe that we can make a nearest neighbor 
grouping jointly for all at, so we can apply Lemma [T2l item 6, to see that both 
denominator and numerator of the second fraction in (|45[) arc continuous. More- 
over, the denominator is > 1, so ind (at) is a continuous function of t. On the 
other hand, numerator and denominator of ind a have to divide every cell dimen- 
sion, so there is only a finite choice of possible values, given the neighborhood 
structure. Hence ind (at) is constant. □ 
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In order to state the remaining Lemma, let us consider any automorphism a 
of a tensor product Al <£> Ar of two finite dimensional matrix algebras. We can 
take A x = B(H X ), for some (^-dimensional Hilbert spaces H x ). In that case, the 
automorphism a is implemented by conjugation with a unitary: a(x) = UxU* , 
for some U E U(Hl <£> Hr)- We will express the fraction appearing in (|45p in 
terms of U. 

In the following calculation, it turns out to be convenient to introduce or- 
thonormal bases in the Hilbert spaces involved. Denote these bases by {\i)} C %l 
and {\a}} C T-Lr respectively. These bases allow us to define the notion of a par- 
tial transpose U r of U. We set 

(ia\(U r )\jb) = (ja\U\ib). 

While the partial transpose depends on the basis used to define it, one may easily 
convince oneself that the expression appearing in next lemma is independent of 
that choice. 

Lemma 14 When a(x) = UxU* is an automorphism of the tensor product 
Al <8> Ar of finite dimensional matrix algebras, we have 

r,(a(A L ),AR) 2 = ^-Tr((U r U r *j 2 ) (46) 

Moreover, this expression is invariant under the substitution Al <-> Ar. 

Proof. We introduce matrix units = \i)(j\ as above and use that the Vdeij 
form an orthonormal basis of Al- Thus the projection Pl on T-L T corresponding 
to this algebra is Pl = dL^2ij l e «')( e yl- Of course, the transformed projection 
Pl for o(Al) is obtained by substituting \a{dj)) for ley). Similarly, we set 
fab = \a){b\ and, accordingly, P R = d R J^ a b \ fab){fab\- Then the left hand side of 
(|46|) becomes 

r,(a(A L ),AR.) 2 = Tt(P l Pr) = d L d R ^ KUHeij)^ 2 , 

ijab 

with the scalar products 

(/a&Key)) 



Altogether we get Eq. (|4T)| . 

The second claim becomes a simple corollary: we must verify that the right 
hand side is unchanged by the substitution U <-> U* . Since the adjoint operation 
commutes with the partial transpose, this amounts to a cyclic rearrangement 
under the trace. □ 



y^^r(e k kfbaUejjfccU*) 

kc 

l —Y,(HU\ic)(jc\U*\kb) 



dhd R 



kc 



)-^2(ia\U r \kc)(kc\U r *\jb) 



dLd R 



dLdR 



{zap 1 U l *\jb) 



(47) 
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8. Outlook 

Two directions of generalization of the theory presented here are especially sug- 
gestive and arc, in fact, the subject of ongoing work in our group. We briefly 
comment on the prospects. 

8.1. Approximate causality. In many situations in physics causality is only ap- 
proximately satisfied, e.g., as a bound ||C4;y|| < f(\x — y\) for some function / 
going to zero at infinity. For example, the unitary groups generated by finite 
range Hamiltonians would satisfy this, but are never strictly causal in the sense 
required here. For the key Lemma. 2] powerful generalizations exist [3]. From 
these it is clear that the index of approximately causal unitarics is definable, 
integer valued, satisfies a product formula, and is zero for the unitaries arising 
from Hamiltonian subgroups. The part of the theory which is less clear is related 
to the converses, namely the construction of strictly causal unitary walks with 
the same index, approximating a given approximately causal unitary. 

In the case of cellular automata the right notions of approximate causality 
are not clear. Ideally, one would only demand that a be an automorphism of the 
quasi-local algebra [5] . By definition, this means that the image of any localized 
element can be approximated in norm by localized ones. The idea of support 
algebras (which largely relies in its present form on picking a finite basis) is 
certainly too simplistic, and in any case it is unlikely that such algebras will 
always turn out to be finite dimensional matrix algebras. One may speculate, 
whether the index should take on also irrational values, but this seems unlikely, 
because of its dependence on the cell structure: for chains of homogeneous cell 
dimension 2, the index is always a power of 2, and not a dense set of rationals. 

8.2. Higher lattice dimension. As we have shown, in lattice dimension 1 three 
possible classifications of walks and automata coincide: (1) the classification 
by locally computable invariants, (2) the classification modulo locally imple- 
mentable, unitaries, and (3) the classification up to homotopy. It is very unlikely 
that these three coincide in higher dimension as well. In the one-dimensional 
theory we could allow the local systems to grow, but also the localization re- 
gions. In the higher dimensional case we will use a translation invariant metric to 
bound the neighborhood sizes of a "local" system. Generalizations can be built 
on coarse geometry [Tl] . 

Locally computable invariants are probably trivial. For example, an arbi- 
trarily large patch of the shift automorphism can be connected to the identity 
outside a finite enlargement of the patch. In this sense the shifts have the same 
invariants as the identity. To get a more interesting theory, one should take other 
regions for the definition of "locally" computable, e.g., computability on cones 
[B] , or computability outside of any arbitrarily large region. 

The classification modulo local implementability is especially interesting from 
the physical point of view, but it might turn out to be rather wild. For example, 
consider some self-intersection free path in the lattice, which comes from infinity 
and goes to infinity. Fix walks/cellular automata, which are equal to the iden- 
tity off this path and allow any of the onc-dimcnsional systems along this path. 
As long as we fix the path, we can apply the one dimensional theory. For two 
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paths, which keep a finite distance from each other, it is easy to envisage local 
swap-type unitaries taking one path to the other, which would bring the corre- 
sponding path-related indices under the same roof. However, if the paths move 
away from each other, there will be no such local operation connecting them, so 
systems with non-trivial indices along these paths fall into different equivalence 
classes modulo local implemcntability. However, the equivalence classes of paths 
modulo "keeping a finite distance" are a rather unmanageably large set. A useful 
classification cannot be expected. Incidentally, the same class of examples shows 
that the "invariants computable outside any finite region" will give a wild set. 

For the homotopy classification of walks there is already a theory, based on 
the K-theory of C*-algebras and its connection with coarse geometry. Indeed, 
the Ki group of a C*-algebra just classifies the connected components of its 
automorphism group. This theory will most naturally apply to approximately 
causal walks, since strict causality cannot even be stated simply in terms of the 
C*-algebra of the whole system. The connection with coarse geometry is being 
explored, for example, by Ralf Meyer and his group in Gottingcn. Surprisingly 
(to us) it turns out that using Bott periodicity one can see that the K\ group 
of approximately causal walks alternates between Z in odd dimension and in 
even dimension. Unfortunately, this theory does not apply readily to the case of 
cellular automata. 



8.3. Higher dimensional translation invariant systems. In order to tame the 
wildncss indicated in the previous subsection, one can restrict attention to trans- 
lationally invariant systems. Immediately, the index of walks gets an obvious 
definition. The Fourier transform U of the walk is now a Laurent polynomial in 
the variables exp(ipk), where (pi, . . . ,p s ) is the momentum vector. Then det U 
is an invcrtiblc polynomial and we again conclude dct U(p) = cxp(i^2 k nkPk) 
for some integers n^. The lattice vector (ni, . . . ,n s ) can be called the index, in 
direct generalization of the one-dimensional case. 

For cellular automata, which are products of partial shifts and local block 
unitaries, we can just define the index via the shift content contained in such a 
representation, obtaining some vector with rational components. Since it is not 
known, whether any QCA is of this form, this is very unsatisfactory. As a step 
in the right direction, one can define an index by reduction to the ID case, with- 
out using any special decomposition: Suppose we choose vectors a\, . . . , a„_i in 
an n-dimensional lattice, and we identify sites differing by integer multiples of 
these vectors. Then if the are large enough with respect to the interaction 
length of the automaton, this gives a well-defined cellular automaton evolution 
a[ai, . . . ,a„_i] on the quotient lattice, which now has only one unbounded di- 
rection, so we can assign an index to it. Then we call a vector q € R™ the index 
of a, if 

inda[ai, . . . ,a n _i] = det[<j,ai, . . . ,a n -i], (48) 

where on the right hand side the square bracket denotes the matrix with the 
specified column vectors. It can be verified easily that this gives the result indi- 
cated before, for any cellular automaton, which has a decomposition into partial 
shifts and partitioned unitaries. But in general it is not even clear that the left 
hand side must depend linearly on each afe. 
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